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Abstract 

Thermodynamical properties of nuclear matter at sub-saturation 
densities were investigated using a simple van der Waals-like equation 
of state with an additional term representing the symmetry energy. 
First-order isospin-asymmetric liquid-gas phase transition appears re- 
stricted to isolated isospin-asymmetric systems while the symmetric 
systems will undergo fragmentation decay resembling the second-order 
phase transition. The density dependence of the symmetry energy 
scaling with the Fermi energy satisfactorily describes the symmetry en- 
ergy at sub-saturation nuclear densities. The deconfinement-confine- 
ment phase transition from the quark-gluon plasma to the confined 
quark matter appears in the isolated systems continuous in energy 
density while discontinuous in quark density. A transitional state of 
the confined quark matter has a negative pressure and after hadroniza- 
tion an explosion scenario can take place which can offer explanation 
for the HBT puzzle as a signature of the phase transition. 

Introduction 

The knowledge of the phase diagram of nuclear matter is one of the princi- 
pal open questions in modern nuclear physics with far reaching cosmological 
consequences. Detailed investigations have been carried out in the recent 
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years in particle-nucleus and nucleus-nucleus collisions in a wide range of 
projectile energies. 

The process of multifragmentation was investigated at intermediate and 
high energies in order to study the properties of the expected liquid-gas 
phase transition at sub-saturation nuclear densities. For instance, using the 
calorimetry of the hot quasi-projectile nuclei formed in the damped nucleus- 
nucleus collisions (see Ref. [TJ and Ref. [2] for a review of related methods), 
thermodynamical observables of the fragmenting system such as temperature 
and chemical potential were extracted and a set of correlated signals of the 
isospin-asymmetric liquid-gas phase transition was observed. On the other 
hand, various properties of the observed fragments can be interpreted as 
signals of criticality [31 SI EJ , thus indicating an underlying second order 
phase transition. The phase diagram at sub-saturation densities appears to 
be rather complex and detailed investigations are needed for clarification. 

The recent experimental studies of the nucleus-nucleus collisions at ultra- 
relativistic energies provide vast amount of experimental observations which 
can be interpreted as signals of production and decomposition of the new 
state of matter with deconfined quarks, called quark-gluon plasma (QGP) 
[6]. Such matter exhibits many unique properties such as nearly perfect fluid 
behavior [7] and strong medium modification of high-energy particles with a 
typical Mach cone behavior [8]. One of the phenomena not yet understood 
there is the so-called "HBT puzzle" [9], an unexpected behavior of the source 
radii extracted from particle-particle correlations, which appear to be almost 
isotropic, in contrast with intense elliptic flow observed in such collisions. 
It was suggested that spatial inhomogeneity (typically in the form of QGP 
droplets) can play a crucial role in its explanation [10]. The mechanism of 
formation of such inhomogeneity is one of the interesting open questions. 

In the present work, the mechanism of a phase transition is investigated 
as a phenomenon reflecting the properties of the underlying interaction. A 
modified van der Waals equation of state is used for the nuclear medium 
at sub-saturation densities and the role of isospin asymmetry will be in- 
vestigated within the limitation following from the evolution of an isolated 
system. An analogous treatment is used for schematic description of the 
system undergoing deconfinement-confinement phase transition. 
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Figure 1: Density dependence of a) pressure and b) internal energy, given 
by the nuclear equation of state from Eq. HJ The lines represent isotherms 
with temperatures increasing with the step 2.5 MeV, starting at zero. 

Isospin asymmetric liquid-gas phase transition 
in the nuclear matter 

The liquid-gas phase transition in the nuclear matter is supposed to occur 
after nuclear matter expands and medium becomes unstable toward density 
fluctuations. The process as initially suggested in Ref. [TTJ was supposed to 
be analogous to liquid-gas phase transitions in real gases, which is typically 
described by the van der Waals equation of state (EoS). However, unlike the 
van der Waals gas, nuclear matter is a two-component system and it thus 
possesses an additional degree of freedom related to proton and neutron con- 
centrations. Additional energy term representing the symmetry energy has 
to be taken into account. Under assumption that the symmetry energy de- 
creases with density, it can be favorable for the expanding system to store the 
excess of neutrons (protons) into a dilute phase and thus isospin-asymmetric 
liquid-gas phase transition can occur [12J. However, the expanding nuclear 
system is isolated and no heat can be transferred into or out of the system 
and a typical picture of phase transition under conditions of phase equilib- 
rium, as represented by the well-known Maxwell construction, is not relevant 
due to necessity to transfer latent heat into or out of the system. 

In order to study the mechanism of isospin asymmetric liquid-gas phase 
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transition in the isolated nuclear system, we adopt a simple EoS for the 
symmetric nuclear matter, closely reminding the van der Waals EoS, where 
the repulsive short-range cubic interaction term is introduced as a substitute 
for the minimum volume 



The parameter values a = —263.7 MeVfm 3 i and b = 1674.94 MeVfm 6 
were chosen to fix zero pressure at saturation density taken as po — 0.16 fm -3 
and preserve plausible values for both critical temperature (T c ~ 14 MeV) 
and ground state energy {U gs ~ —21 MeV). The equation of state thus 
can be considered as a reasonable approximation to the results obtained 
using the microscopic theory of nuclear matter. The resulting pressure and 
internal energy isotherms are shown in Fig. [TJ The internal energy (obtained 
assuming the heat capacity cy = 3/2) exhibits a quadratic dependence on 
density with a minimum at saturation density, what is the expected result 
for the nuclear ground state. 

Since the system is thermally isolated, the expansion of the system is an 
adiabatic process. The entropy does not change during system evolution and 
assuming that the heat capacity cy is constant, the following expression for 
temperature as a function of density 



can be derived, where To is the temperature projected to saturation den- 
sity po- The resulting isoentropic density dependences of the pressure and 
internal energy are shown in Fig. [2^, b. One can identify spinodal region 
where pressure grows while density falls, however, the limiting temperature 
appears much smaller than the one extracted from Fig. [TJ The resulting 
phase diagram constructed using the pressure and temperature at the spin- 
odal contour is shown in Fig. [2b. The critical temperature appears to only 
slightly exceed 8 MeV. Such observation is in agreement with experimen- 
tal results [13]. Two caloric curves are shown in Fig. [2H. The dashed line 
represents the standard caloric curve obtained using the Fermi-gas formula 



where the level density parameter is chosen as a = ^4/(10 MeV). Such 
dependence is widely used in various evaporation models and represents es- 



p = pT + ap 2 + bp 3 



(1) 





E* = aT 2 , 



(3) 
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sentially an evolution at ground state density. It is worthwhile to note that 
the temperature dependence in Eq. [3] is equivalent to the heat capacity cy 
linearly increasing with temperature. Such dependence is in agreement with 
theory of the Fermi gas at low excitation energies, however, its validity should 
be restricted only to the domain where cy does not exceed the ideal gas value 
3/2 (below temperature 7.5 MeV), since cy saturates at this value and thus 
assumption of its further growth is not correct. Solid line in Fig. [2U repre- 
sents the caloric curve obtained by combining the temperature at which the 
system enters spinodal region with the estimate of initial excitation energy 
taking into account also expansion of the system during initial heating. The 
heat capacity c p was obtained from 



c p = °v + m ; — , ^T3. ' ( 4 ) 



T 

T + 2a~p~+ 3bp 2 

where cy was determined using the formula [3J The system expanded 
according to 




= T + lap + 3bp 2 ' ^ 

which is closely related to c p . At each step of the procedure the tem- 
perature was incremented by a small value, corresponding increase of the 
internal (excitation) energy was proportional to c p and the density was mod- 
ified according to the calculated volume expansion coefficient. Finally, new 
value of the pressure was calculated from the equation of state using new 
values of density and temperature. The resulting evolution of the system is 
shown in Fig. [2^ as a nearly horizontal line starting from the ground state. 
The pressure of the system appears to grow very slowly while the system 
expands. The resulting caloric curve is similar but much flatter than the one 
obtained using the Fermi gas formula. For real nuclei one can assume that 
the system will follow the dashed line until it has enough energy to pass the 
fragmentation barrier and expand toward spinodal region, thus approaching 
the solid line via a short plateau. Such a short plateau was observed in the 
experiment [T], where it correlated with the onset of thermal equilibrium 
with the nucleonic gas, as represented by agreement of the extracted double 
isotope-ratio temperature with the kinematic temperature of protons. 

Furthermore, it is interesting to note that while the fermionic system ex- 
pands its Fermi energy exhibits the p 2 / 3 density dependence, which is iden- 
tical to density dependence of temperature during adiabatic expansion when 
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one assumes the saturation value cy = 3/2. Thus the ratio of the Fermi 
energy to temperature, which is a parameter determining the occupation 
numbers, will remain the same during adiabatic expansion. The ideal gas ap- 
proximation cy = 3/2 appears justified for the isoentropic expansion where 
such value was reached during heating, in the Fig. [2^ such situation applies 
to pressure curves with initial temperatures above T = 9 MeV. 




p (MeV/fm 3 ) E* (A MeV) 



Figure 2: Density dependence of a) pressure and b) internal energy for 
constant entropy. The lines represent adiabatically expanding systems with 
the temperature at saturation density T increasing with the step 2.5 MeV, 
starting at zero. Curve in the panel a) represents the spinodal contour and 
the horizontal line heating process of the nucleus. In the panel c) is the p-T 
diagram for the spinodal contour and in the panel d) are the caloric curves 
for the spinodal contour (solid) and for the Fermi gas formula (dashed). 

Nuclear matter is a two-component system and thus an energy term rep- 
resenting the symmetry energy has to be introduced in order to reflect the 
difference of proton and neutron concentrations 
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where a a is the symmetry energy coefficient from the Weiszacker formula 
and a p 2//3 density dependence is assumed from the behavior of the Fermi en- 
ergy. Due to density dependence of the symmetry energy it can be favorable 
for the system to store the excess of neutrons ( protons ) into a dilute phase 
and thus the isospin-asymmetric liquid-gas phase transition can take place. 
The equation of state will be modified to the form 

rp | 2 . , 3 , 2a aP0 ( P N5/3 ( N ~ Z V / 7 N 

P = P T + ap +b P {—) , (7) 

which allows us to investigate the behavior of the system as a function of 
density and isospin asymmetry / = (N — Z)/A. The additional term in the 
equation of state will not modify the temperature dependence in Eq. [2J 




Figure 3: a) Pressure contour plot, corresponding to the nuclear equation 
of state given in Eq. [7J The adiabatic expansion with the temperature at 
saturation density T = lOMeV is assumed. Thick curve represents the 
spinodal contour. In panel b) is shown a detailed view onto low density 
region. 

The pressure contours determined by Eq. [7J using the value a a = 25 MeV 
are shown in Fig. [3] as a contour plot. The resulting spinodal region is 
highlighted and its shape is similar to microscopic calculations [H]. The 
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asymmetric systems will enter the spinodal region at lower densities, on the 
other hand the isospin-asymmetric gas remains stable at much higher densi- 
ties than the symmetric one. 




Figure 4: The resulting values of a) proton concentrations x v (solid line), Xi 
(dashed line) and b) vapor fraction a vap for an adiabatically expanding system 
with T = 10 MeV and isospin asymmetry of the system (N — Z)/A = 0.3. 

Metastable state inside the spinodal contour ((dp/dV)s > 0) is unsta- 
ble toward spontaneous compression to stable liquid state at the spinodal 
contour, but such process would create heat. Instead of heat, the released 
internal energy can be deposited into expansion of vapor into free volume, 
thus spending released energy into overcoming attractive nuclear force. When 
the whole system is isolated (5Q tot = 0) the process of spinodal decomposi- 
tion will be equivalent to fast transfer of particles from the metastable state 
at its constant pressure to the stable liquid and vapor states at their corre- 
sponding pressures and total enthalpy of the system will not change during 
phase transition and the relation 

U ms + PmsVtot = Uu q + PliqVuq + U vap + P V apV vap (8) 

will be preserved. 

During fast transition the total volume remains constant (dV to t = 0) and 
for a given combination of the stable liquid and gas states located at the 
spinodal contour the partial volumes of liquid and vapor are 
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Figure 5: The resulting values of a) proton concentrations x v (solid line), x\ 
(dashed line) and b) vapor fraction a vap for an adiabatically expanding system 
with T = 10 MeV and isospin asymmetry of the system (N — Z)/A = 0.03. 
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where 



Hot -< *-vap 

and nucleonic fraction of the vapor is 

A, 
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(10) 



(11) 



Hot Pms Pliq Pvap 

for which the charge conserving proton concentration (x = Z/A) is de- 



termined using proton concentration of a given liquid state xi 



x 7 



Xn\p_ 



&vap)%l 



(12) 



l vap 



which should agree with the proton concentration x v of the vapor state 
at the spinodal contour used to determine a vap . Such a self-consistent stable 
liquid and vapor states located at the spinodal contour can be determined 
numerically. Total enthalpy conservation relation can be written as 
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H m s Hi H v h ms (Lvaphv (1 Q"vap)hl , 

where h ms , hi and h v are values of the enthalpy per particle. 
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Figure 6: The resulting values of a) proton concentrations x v (solid line), xi 
(dashed line) and b) vapor fraction a vap for an adiabatically expanding system 
with Tq = 10 MeV and isospin asymmetry of the system (N — Z)/A = 0.003. 

The resulting values of x v , xi and a vap for an adiabatically expanding 
system with T = 10 MeV are shown in Figs. HJ El and M for the isospin 
asymmetries of the system (N — Z)/A = 0.3, 0.03 and 0.003, respectively. 
It can be seen that for the asymmetric system with (N — Z)/A = 0.3 there 
occurs a well pronounced isospin asymmetric liquid-gas phase transition with 
a sizable part of the system being converted into very isospin asymmetric 
vapor while the rest of the system becomes still more and more symmetric 
liquid. These trends become more prominent as the density of initial meta- 
stable state decreases. When the system becomes more symmetric (see Figs. 
|5] and [6]) the vapor fraction drops very quickly while the vapor preserves its 
isospin asymmetry. Also the region of densities where isospin asymmetric 
liquid-gas phase transition can occur shrinks quickly. From the Figs. IH 
and [6] one can conclude that observable isospin asymmetric liquid-gas phase 
transition will occur for very asymmetric systems such as these observed in 
the experiment [1] while for symmetric systems (and for the remnants of 
liquid phase in asymmetric systems) another mechanism will play role. In 
order to reach spinodal region the system must pass the multifragmentation 
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barrier and thus these systems will be unstable toward sudden decomposition 
into multiple fragments. Such a scenario can provide explanation for the 
somewhat confusing situation where signatures of both first and second order 
phase transitions can be extracted from observed data (even from the same 
set of data). The first-order behavior can be attributed to the distillation of 
isospin-asymmetric gas while the second-order signatures can be attributed 
to subsequent fragmentation (percolation) of the remaining symmetric liquid. 

It is also interesting that the conclusions observed here are in good agree- 
ment with the assumptions used by the Statistical Model of Multifragmen- 
tation (SMM) [15] where the fragment partition is considered at equilibrium 
with nucleonic gas. Such a scenario is supported also by experimental data 
PP. Other interesting point is that the p 2 ^ 3 density dependence of the sym- 
metry energy used in the present work appears as satisfactorily describing 
the properties of asymmetric nuclear matter at sub-saturation densities. 

Deconfinement-confinement phase transition in 
the quark matter 

The ultrarelativistic nucleus-nucleus collisions result in the very high energy 
densities where the quark matter can exist in thermally equilibrated decon- 
fined state known as quark-gluon plasma. During the subsequent expansion 
energy density decreases until the critical value around 1 GeVfm~ 3 where, 
according to the predictions of the quantum chromodynamical calculations 
performed on the lattice [16] the matter will revert to the confined state, 
where the interaction is dominated by a divergent confinement potential, 
linearly increasing with the distance. The potential energy, stored in the 
one-dimensional field tubes, called strings, will be subsequently released by 
string fragmentation in the form of hadronic gas. In order to investigate the 
phase diagram, one can use known equations of state of initial and final states 
of the quark-gluonic plasma and the hadronic gas and try to interpolate the 
equation of state in the spinodal region, as is done in Ref. [T7]. However, 
since the process of hadronization typically proceeds via fragmentation of 
elongated strings it is also interesting to explore thermodynamical properties 
of the possible intermediate state of matter, where the particles of the quark 
gas interact via confinement potential. Such state of matter is referred here 
as the confined quark matter. The properties of intermediate states can be 
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important for the mechanism of phase transition, as was demonstrated in the 
case of liquid-gas phase transition at sub-saturation densities. Here we will 
explore thermodynamical properties of the confined quark matter in analo- 
gous way as thermodynamical properties of the asymmetric nuclear matter 
in previous section. 

In order to investigate properties of the confined quark matter a confine- 
ment potential is defined as 

Vconfir) = nr , (14) 

where k is the string tension with a typical value of 1 GeV/fm, as used 
in the string fragmentation models [32]. In order to explore properties of 
the confined quark matter, a 3-dimensional cubic lattice with a quark-quark 
distance L can be assumed, where each quark interacts with six nearest sites. 
Then the component of the internal energy, corresponding to confinement 
potential, can be determined as 

U conf = ^N q (6nL) = SNgnp- 1 / 3 , (15) 

where the number of quarks N q is large enough to neglect surface effects. 
However, an assumption that quark on the lattice interacts with all neighbors 
must not necessarily be satisfied, since strings are essentially one-dimensional 
objects. One can modify above formula by assuming that each quark on the 
lattice interacts with w neighboring quarks and one obtains the formula 

U conf = \N q {wKL) = N^p- 1 / 3 . (16) 

The pressure component corresponding to confinement potential can be 
obtained as 

A _ dUconf W 2/3 1 Uconf n? , 

thus leading to negative pressure. In order to obtain the equation of state 
of the confined quark matter, one has to introduce the pressure of the quark 
gas which in the ultrarelativistic limit can be written as 
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where g q is a degeneration factor and r\ is the Dirichlet Eta function. In 
order to write equation of state with the confinement potential one needs to 
find relation between energy density and quark density That is possible for 
zero chemical potential (with g q accounting for equal numbers of quarks and 
anti-quarks) 



P 



N 
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h 3 



t 3 v (3) 



and after comparing the equations ([15]) and (|T9|) one arrives to 

A*) 



p = pT 



for the equation of state. 
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Figure 7: Dependence of the pressure on a) energy density and b) quark 
density for w = 2. Thick line represents quark-gluon plasma (Eqs. (Il8j) 
and (1T9|1 with g q = 24 and g g = 16), solid section corresponds to the region 
above phase transition. Thin solid line shows pressure of the confined quark 
matter, dashed line shows the quark gas contribution (Eqs. (|T8|) and (|T9l 
with g q = 24), dash-dotted line shows the confinement potential contribution 
(Eq. (ED). 

The resulting situation is shown in Fig. [7] for w — 2. At energy densities 
below critical value 1 GeV/fm 3 and corresponding quark densities the con- 
fined quark matter exhibits negative pressure and also violates the mechan- 
ical stability criterion. Such a medium is indeed unstable against collapse 
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what on the other hand is a direct consequence of the interaction rather 
than an interplay of short- and long-range interactions, as in the liquid-gas 
phase transition. Since confined quark matter is primarily unstable toward 
hadronization it can be assumed that the survival time against hadronization 
is much shorter than collapse time and the medium can be considered as a 
quasi-stable transitional state. 

The detailed mechanism of the deconfinement-confinement phase transi- 
tion is beyond the scope of this work, the lattice QCD calculations [16] sug- 
gest a rapid crossover for the systems with small chemical potential produced 
at RHIC. As in the liquid-gas phase transition at sub-saturation densities, 
one can explore the global properties of the interaction and limitations on the 
phase transition resulting from the isolated system and, as a consequence, an 
isoenthalpic transition into confined quark matter can be considered. In the 
present case the hot quark-gluon plasma with high thermal pressure trans- 
forms into confined quark matter with much lower (negative) pressure. The 
hadronization phase transition then follows as a consequence of the fragmen- 
tation of the elongated strings. 

Enthalpy of the quark-gluon plasma with u and d quarks and antiquarks 
can be expressed as 

where g q = 24, g g = 16 and ( is the Riemann Zeta function. The en- 
thalpy balance of the transition from quark-gluon plasma at temperature T 
to confined quark matter at density p and corresponding temperature T can 
be written assuming ultrarelativistic fermionic gas also in the confined phase, 
with the additional component corresponding to confinement potential, ob- 
tained after combining the equations [16] and [T7] 

*N,m + «^)T„ = 4N,'Mt + 2N "^V (22) 
7?(3) g q rj{3) 77(3) 6 

and using the relation [19] one arrives to equation 

4(^§ + ^)r„ = M . h + ^ (23) 

which after assuming numerical values T = 170 MeV (as suggested in 
the Refs. [6], PH), k = 1000 MeV/fm p], g g = 24 and g q = 16 leads to 
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Figure 8: Solutions of the Eq. [231 Dependence of the a) parameter w, b) 
energy density, c) pressure and d) temperature of the confined quark matter 
on quark density. Solid lines represent the low-density branch, dashed lines 
show the high-density branch. 

a numerical equation which has solutions for w < 1.862. The solutions of 
Eq. (I2"3"j) are shown in the panel Fig. [S^ as a dependence of the parameter w 
on quark density. A low-density branch (solid line) and high-density branch 
(dashed line) can be identified, separated by the point at p — 0.921 fm -3 
corresponding to maximum value of the parameter w = 1.862. The resulting 
dependence of energy density on quark density is shown in the panel Fig. [8b 
and one can see that the solution for w = 1.862 corresponds approximately 
to initial energy density of quark-gluon plasma at To = 170 MeV. The 
low-density branch corresponds to lower energy densities and thus represents 
possible states to which the system can transform since it corresponds to both 
quark and energy densities below the phase transition. The high-density 
branch represents the states with higher energy densities than the critical 
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density and thus can be considered unphysical since it is not obvious that 
the confinement potential is applicable to such states. From the Eq. 0231) 
thus implies that the phase transition can be in principle continuous in the 
energy density, but discontinuous in the quark density. After reaching the 
critical energy density, further expansion will be achieved by transforming 
an increasing amount of the quark-gluon plasma into confined quark matter 
with the maximum of about two thirds of quark density, temperature below 
150 MeV and negative pressure (see Fig. [Hb). However, once the strings will 
fragment and a hadronic gas will be formed, potential energy will be released, 
pressure will increase dramatically and one can assume an explosion scenario. 
One can imagine a sudden expansion of the bubble of hadronic gas to the 
radius reaching the dimensions of the whole system. Such a scenario can 
offer a natural explanation to the "HBT puzzle" where the extracted radii 
are essentially equal for longitudinal, sideward and outward direction and 
suggest an isotropic source with a short emission time. HBT puzzle can thus 
be interpreted as a signature of the phase transition. It is also interesting 
to note formal similarity with the cosmologic inflation, a scenario leading 
to isotropic causality-connected universe, which is supposed to start from 
the patch of matter in the state of false vacuum at negative pressure. In 
the present case the confined quark matter could represent the false vacuum 
while the hadronic gas could be considered as the true vacuum. Thus, the 
ultrarelativistic nucleus-nucleus collisions may exhibit, apart from the "Little 
Bang" scenario [TTJ] and possibility of the existence of mini black holes [20J, 
another cosmo logical analogy in the form of "mini-inflation", as a possible 
explanation of the HBT puzzle. 

It is also interesting to compare the results for the quark matter to the 
liquid-gas phase transition at sub-saturation density. Both systems will 
evolve spatial inhomogeneity with one phase at positive pressure and the 
other one at negative pressure, but different asymptotic behavior of the po- 
tential, one diverging at small distances while the other one at large dis- 
tances, leads to reversion of some properties. While in the liquid-gas phase 
transition the dense liquid phase transforms into dilute gaseous phase, in the 
deconfinement-confinement phase transition the dense gaseous phase trans- 
forms into dilute liquid phase, which subsequently transforms into hadronic 
gas. The liquid consistence of the confined quark matter can be judged 
e.g. from its possible one-dimensional structure, as suggested by existing 
isoenthalpic solutions. In both cases the liquid phase is unstable toward 
subsequent decay by fragmentation mechanism. It is also interesting to in- 
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vestigate whether introduction of strangeness into the model description will 
lead to asymmetric phases. One can assume that for the strange quarks the 
string tension will increase and thus it will be preferable for strange particles 
to remain in the deconfmed phase. The strangeness of the deconfined phase 
will increase and the strange hadrons will be formed in the latest stage of 
hadronization. 

Summary and conclusions 

Thermo dynamical properties of the nuclear matter at sub-saturation den- 
sities were investigated using a simple van der Waals-like equation of state 
with an additional term representing symmetry energy. For the isolated sys- 
tem an enthalpy conservation rule was introduced resulting in significant 
limitations for the isospin-asymmetric liquid-gas phase transition, which is a 
sizable effect for very isospin-asymmetric systems but becomes negligible for 
symmetric systems. First-order phase transition thus appears restricted to 
asymmetric systems while for symmetric systems and symmetric remnants 
of the first-order phase transition the multifragmentation decay can exhibit 
similarity with second-order phase transition. The density dependence of the 
symmetry energy scaling with density dependence of the Fermi energy satis- 
factorily describes the symmetry energy at sub-saturation nuclear densities. 

The limitations resulting from the isolated system were investigated also 
for the deconfinement- confinement phase transition from the quark-gluon 
plasma to quark matter with the confinement potential. It appears that the 
transition can be continuous in energy density while discontinuous in quark 
density. A transitional state of confined quark matter has a negative pres- 
sure and explosive scenario can take place during hadronization. Resulting 
scenario can offer an explanation for the HBT puzzle as a signature of the 
phase transition. 

This work was supported through grant of Slovak Scientific Grant Agency 
VEGA-2/5098/25. 
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